Current-induced motion of a domain wall in magnetic nanowires 
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The dynamics of current-induced motion of a magnetic domain wall in a quasi-one-dimensional 
ferromagnet with both easy-axis and easy-plane anisotropy, is studied. We pay a special attention 
to the case of a sharp domain wall, and calculate the spin torque created by the electric current. 
The torque has two components, one of which is acting as a driving force for the motion of the 
domain wall while the other distorts its shape, forcing thus the magnetic moments to deviate from 
the easy plane. 

PACS numbers: 75.60.Ch,75.70.Cn,75.75.+a 
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I. INTRODUCTION 

The increased interest in the dynamics of domain walls 
in magnetic nanostructures is rooted mostly in perspec- 
tive applications in novel spintronic devices. Several re- 
cent experiments demonstrated that the motion of do- 
main walls can be effectively controlled by means of an 
external magnetic field or an electric current iS* 3 ^ 4 ^ 

The theory of domain wall motion has been worked 
out long time agOi&L& Within this theory, the magnetic 
dynamics is usually described in the framework of a clas- 
sical ferromagnet assuming separation of the magnetic 
and electronic degrees of freedom. In this way three- or 
two-dimensional ferromagnets are treated. 

Present studies of the domain wall dynamics are de- 
voted to the issue of the current-induced magnetic mo- 
tion in nanowires and nanoconstrictions. In the presence 
of an electric current, the domain wall can move due to a 
spin torque exerted on the magnetic system by the spin- 
polarized electron gas. In addition, (linear) momentum 
can be transferred directly to the domain wall upon scat- 
tering of the charge carriers. The problem of calculating 
the spin torque is of a prime importance for the theory of 
the domain-wall motion, however only few recent works 
address this point specificallyi^iiMiiiSiii 

The domain-wall motion can be described by Landau- 
Lifshitz equations, which have a well-known static so- 
lution for the domain wall. Finding the corresponding 
dynamical solution is however a nontrivial task, particu- 
larly in the presence of an external force. The standard 
way is to describe the domain wall dynamics using an ap- 
proximate scheme, which is based on physically reason- 
able arguments, but a strict mathematical justification 
is missing. The simplest approach is to use an approxi- 
mate solution that assumes the moving wall to have ex- 
actly the same shape as deduced from the static solution. 
The rigorous conditions for the range of validity of this 



approximation are however unclear. To obtain reliable 
results numerical simulations are needed^ 4 - 

Most recently, new approaches have been put forward 
which address&i^ the calculation of the torque as well as 
the solution of the dynamical equations of motion of the 
domain wall. Some aspects of these works are developed 
further and partially revised in this paper. 

Our treatment is concentrated on the spin torque and 
on the wall dynamics in a magnetic nanowire with a 
domain wall which is sharp on the length scale set by 
the wave length of the relevant charge carriers. The ap- 
proach is thus more appropriate for the case of magnetic 
semiconductors with a small Fermi momentum of carriers 
(electrons or holes) i& 



II. CURRENT-INDUCED SPIN TORQUE 

Here we consider the spin torque transferred by the 
electric current from the spin polarized electron system 
to the domain wall. The main goal of our calculation is 
the demonstration of the presence of two components of 
the torque, which tend to rotate the magnetic moment 
in different directions. 

We adopt a one-dimensional model for the charge carri- 
ers with a point-like interaction between the electron spin 
er and the magnetic moment M(x) located at a point x 
along the wire 



H^t = gv ■ M(x), 



(1) 



where g is the coupling constant. Here the one- 
dimensionality of the electronic system means that we 
consider the electrons within a wire with transversal di- 
mensions smaller than the electron wavelength A, so that 
only the lowest electron subband is relevant. Strictly 
speaking, the coupling of the moment to the electron 
spin depends on the coordinates y, z that characterize 
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the location of the moment within the wire, g(y, z) ~ 
\ipo(y,z)\ , where ipo(y,z) is the wave function of the 
transverse motion of electrons in the lowest subband. For 
simplicity, in the following we neglect this dependence, 
substituting it by an averaged coupling, g(y, z) — > g = 
A~ 4 J g(y,z)dydz, where A is the cross section of the 
wire. 

Considering the scattering of electrons from a magnetic 
moment M(x) we assume that the moment is frozen in 
the point x on the scale of the characteristic times of the 
electron motion. This assumption renders possible calcu- 
lation of the torque as in the case of a static domain wall. 
The calculated torque is then utilized for the investiga- 
tion of the domain wall dynamics. Our assumption relies 
on an adiabatic approximation insofar as we require the 
time scale for the motion of the magnetic subsystem to 
be slow as compared to that of the electrons. 

To calculate the torque in the case of a sharp do- 
main wall, we start from a model describing the one- 
dimensional scattering from a localized moment in non- 
magnetic and magnetic wires. Then we use the results 
obtained for the simplified models, to calculate the torque 
acting locally on the moments within the domain wall. 



A. Single magnetic moment in a nonmagnetic wire 

Let us consider first scattering of electrons in 
a nonmagnetic one-dimensional system (nonmagnetic 
nanowire). The electrons are scattered from a single 
frozen magnetic moment Mo situated at the point x = 0, 
i.e., M(a;) = M.q5(x). Here we denote the coordinate 
along the wire as x. In the absence of the spin-orbit inter- 
action, it is more convenient to use a different coordinate 
system (x',y',z') for the spin space. We calculate the 
total scattering amplitude (beyond the Born approxima- 
tion) of an electron with an arbitrary spin polarization, 
coming from x = — oo, and elastically scattered into the 
states with different spin polarizations. 

Assuming the quantization axis z' along the moment 
M , we can write the spinor wave function of electrons 
as 



ip(x) 



. ikx 



'ikx 




+ e 



-ikx 



ra 
r*b 



x < 0, 



(2) 



x > 0, 



where the coefficients a and b correspond to an arbitrary 
spin polarization in the incident electron wave, r and 
t are, respectively, the reflection and the transmission 
coefficients in the spin up channel 



r = — - 



1 



t 



1 



1 



(3) 



a = gM m/kh 2 , and M is the magnitude of the localized 
magnetic moment. 



Using Eqs. (2) and (3) we can calculate the spin den- 
sity in the wire, associated with the wave function (2) 



Sfa) = ^(x)a fl ip(x). 



(4) 



Sfj, (x) stands for the oscillating spin density created by 
the electrons reflected backwards from the localized mo- 
ment, x < 0. The spin density components are 

l + a 2 + 2a 4 2a 3 

S X >{X < 0) = S x > n ; 1- V 



(1 + a 2 ) 2 

fa 2 a \ 

2 cos (2kx) {s x , + V J . ( 5 ) 



(1 + a 2 ) 2 

v 2 



Syl(X < 0) = Sy 



1 + a 2 + 2a 4 



2a A 



v (1 + a 2 ) 2 x (1 + a 2 ) 2 
+2 cos (2kx) (-s y , -^—^ + s x , Y^tf) < ( 6 ) 



A + 2a 2 , 1 , 2a 2 

S z ,{x < = 8x , — — T - cos (2kx) — — ^ 
\ I + a A 1 + a z 



sin (2kx) 



2a 



l + a< 



, (7) 



where s M is the unit vector determining the spin polariza- 
tion in the incident wave. Similarly, we obtain for x > 

. . I — a 2 2a 

[X > °> = Sx ' (1 + a 2 ) 2 ~ Sy ' (1 + a 2 ) 2 ' (8) 

n , ~s 1 — a 2 2a 

> 0) - V j^—^ + s x , —^- 2 , (9) 



S z ,{x > 0) = s z , — — - 
1 + a z 



(10) 



As follows from Eqs. (5) to (10), the spin density induced 
by the spin-polarized wave incoming from x = — oo, os- 
cillates with the period n/k at x < 0, and is constant for 
x > 0. 

The spin current is defined as 

3 fa) = ^ [(Vx^W) ^ 1>{x) - ^(x) VMx)] , 

(11) 

(/i = x' , y' } z') and it can be also calculated using Eqs. (2) 
and (3) for x < and x > 0, respectively. Then we find 
that the spin current is constant for x < and x > 0, 
with a jump of x' and y' components at x = 0. 

We can calculate the spin torque acting on the moment 
Mo as the transferred spin current at the point x = 



Ju(S) - fJ+S). 



Using Eqs. (2), (3), (11) and (12) we obtain 



rp Jo 
J-x' — — 

e 



4a 2 2a(l - a 2 ) 

St.' ~ ~ T" Sy' 



1 + a 2 



1 + a 2 



(12) 



(13) 



3 



T , 



JO 



4a 2 



2a (1- a 2 )" 1 

n 



Syl 



1 



(14) 



and T z ' = 0, where e is the electron charge (e < 0), jo is 
the electric current 

(15) 

and w = hk/m is the velocity. Note that we found the 
components of torque (13), (14) in the coordinate system 
related to the moment Mo, so that x' and y' axes are 
perpendicular to the vector Mo, and they are not related 
in any way to the direction of the current jo- 

Using Eqs. (5)-(10) one can show that exactly the same 
values for the torque components (13), (14) are obtained 
by utilizing the relation 



rp 9^0 q ff .s 

= ^— e^x n v &\{Q), 



(16) 



which follows from the equation of motion of the mag- 
netic moment M , where n is the unit vector along M , 
and e^vx is the unit antisymmetric tensor. 

One can also present the result for the torque (13), (14) 
in a form, which is more appropriate for an arbitrary 
coordinate system for the electron spin (not necessarily 
with the axis z along M ) 



= [77 (6^ v - n^n v ) s v + ( e^x n v s\] , (17) 



where 



4a 2 2a (1 - a 2 ) 
'I = ^—2 > C = 7^-^ • (18) 



1 + a 2 



1 + a 2 



As we see from Eqs. (17) and (18), there are two com- 
ponents of the torque - both transverse to the localized 
moment. Apart from this, one of them tends to align 
the moment along the direction of the spin polarization 
of the incoming electrons, whereas the other one is per- 
pendicular to the spin polarization of the incident wave. 
Note that in the Born approximation for the scattering 
amplitude, which is valid for a <C 1, only the second 
term in (17) survives. It rotates the moment Mo to the 
direction perpendicular to vector s (and also to n). 

Thus, the spin torque acting on a single magnetic mo- 
ment can be found as a change of the spin current due 
to scattering from the localized moment. The same re- 
sult is obtained from the calculation of the interaction 
of accumulated spin with the localized moment. In the 
following, to calculate the torque in the domain wall, we 
will use both methods but the second one (coupling to 
the accumulated spin) is more convenient in case of a 
sharp domain wall. 

We assume that in the case of a smooth domain wall, 
it can be more convenient to solve the problem consid- 
ering it as a propagation of a spin-polarized wave with 
subsequent scattering from different magnetic moments. 
Correspondingly, one can calculate the torque as a diver- 
gence of the spin current. 



B. Scattering from a single magnetic moment in a 
magnetic wire 

Now we calculate the torque in the case of a magnetic 
wire with the magnetization M being oriented along the 
axis x for x < (left to the wall) and in the opposite 
direction for x > (right to the wall). Here we as- 
sume the spin coordinate system (x',y',z') to coincide 
with the the (x,y, z) one. Like in the previous prob- 
lem, we introduce an additional frozen magnetic moment 
M = Mo (n x , n y , 0) located at the point x = 0. For 
definiteness, let the vector Mo lie in the x — y plane. 

The relevant Hamiltonian has the form 

h 2 d 2 

H = -^^2 + 9 M(J x s § n T ) + 3 M o n ' trS{x), (I 9 ) 

First, we consider the torque created by a single spin- 
polarized wave (with the spin polarization along the axis 
x labelled as "t") coming from the left. We choose the 
quantization axis along the axis z. Then the wave func- 
tion containing the reflected and transmitted waves of 
opposite polarization is 



v/2 



tr e J- 



>/2 ^ -1 

for x < 0, 



(20) 



V2 



tire T" 



1 



V2 y-i 

for x > 0, 



where fcf^ = [2m(e =F gM)] 1 / 2 /%, and e is the energy. 
Note that the spin-up electrons are the spin-minority 
ones, while the spin-down electrons are the spin-majority 
ones. 

Using the continuity of wave function at x = and the 
condition resulting from the integration of Schrodinger 
equation in the vicinity of domain wall, 



h 2 / dtp 
2m \ dx 



+s 



chp_ 
dx 



+gM (n x cr x + riycry) ip(0) = 0, 

(21) 



we find the transmission coefficients for the spin-up po- 
larized wave 



2fc T (fcj + fcj, - ig n x ) 
(h+hf + g 2 



T/ 



2g Q n y k^ 



(ki+ktf+g 2 



(22) 



(23) 



and the reflection factors r-\ = tf — 1, r-j-/ = t^f, where 
we denote go = 2gmMo/h 2 . 

Using Eqs. (11) and (20), we can calculate the compo- 
nents of the spin current induced by the incoming spin-up 
wave 

w T (1 - |r T | 2 ) + vi |r T /| 2 , x < 0, 
3U X ) = { ( 24 ) 



4 



t f Im [u-f ( e lfc + x — rj e * fe - x ) 



+^4 ( e 



x < 0, 



(25) 



t T /Im — V| £* e zte - x + vi tf e 



-ik—x 



, x > 0, 



states, and the expressions for transmission coefficients 
are 



2ki (fcf + fcj + ig n x ) 
(fc T + fci) 2 +5o 



(31) 



%Re [u T (e ife + a: - r T e"^- 21 ) 

, x < 0, 



i T/ Re [u t if e ife -* + ^ t T e"*-* 



(26) 



x > 0, 



where fc± = fcf ± fcj and v^i = Tik^Jm. 

Hence, the transverse components of the spin currents, 
j^y(x) and j^ z (x), are nonzero for x < and for x > 0. 
As we see from Eqs. (24)-(26), the transverse components 
of the spin current are oscillating functions of x. The 
nonconservation of spin current in the magnetic wire is 
related to indirect magnetic interactions accompanying 
the inhomogeneous distribution of the spin density. In 
the nonmagnetic case, corresponding to the limit of fc_ — > 
0, it reduces to the conservation of spin current at x < 
and x > 0, as we found in the previous section. 

The calculation of the corresponding spin transfer (12) 
in this model gives us the following expressions for the 
components of the torque 

T ]x = 2 V] Retf + (« T + Vl ) (|t T /| 2 - |i T | 2 ) , (27) 



Tiy = -2ti,(vi+v l )Tmti, (28) 



Tf z — 2if / [u-f — (t)-f + f x) Re tf] . (29) 

As in the case of a nonmagnetic wire with a localized 
moment, the torque Tj can be calculated using Eq. (16) 
with the spin density Sj(0), where the index " |" indicates 
the accumulated spin associated with the incident wave of 
spin-up polarization. Using (4) and (20) we can calculate 
the spin density Sf(ar) as in the previous section. The 
spin density oscillates for both x < and x > 0, with 
the periods 2n/k + and 27r/ |fc_| as for the spin current in 
Eqs. (24)-(26). 

With our choice of the coordinate system, when the 
vector Mo lies in the x — y plane, the transverse compo- 
nents of the torque acting on the moment Mq are 



T T j_ = -n y T Vx + n x T-| 



Vy ■ 



(30) 



rotating the moment in x — y plane, and Tj 2 , rotating 
it in out-of-plane direction. The label _L in (30) means 
the projection of torque on the direction perpendicular 
to the moment Mo in the x — y plane. 

Similarly, we can consider the scattering of electron 
with the incident spin polarization opposite to the x-axis 
(labelled as !). The corresponding scattering state has 
the form of Eq. (20) with fcj <-» fcj_ and interchanged spin 



hf 



2g n y k l 
(fct+*l) 2 +flo ' 



(32) 



where the change go — * —go is equivalent to the flip of 
momentum Mo. The components of spin current j s ^ 
have the following form 



-^(i-kiH-filn/l 2 . x <0 > 

-WT|tll 2 +«ll*i/| 2 . ^>o, 



(33) 



-tif Im [vi (e lk+x 



r l e lk - x ) 

, x < 0, 



tjlm vit*e lk - x -v^tie~ lk - x 



(34) 



x > 0, 



t lf Re [v L (e ik + x - ri e- lk - x ) 



-v T (e~ ik + x + rf e lk - x ) 



x < 0, 



(35) 



hf Im 



Vi t\ e ik ~ x +Vi tit 



, x > 0. 



The torque components (27)- (29) can also be calcu- 
lated from Eq.(16), taking into account the spin accumu- 
lation at x = 0, 



Su*(o) 
Su»(o) 
St,i*(o) 



= T 



4fc T 2 a [(fc T +*u) 2 +gg(^-"g)] 
Kfc T +fcl) 2 +.9o 2 ] 2 



pT+fc|) 2 +5o 2 ] 



212 



, (36) 
(37) 
(38) 



The above formulas will be used later to calculate the 
torque exerted on a domain wall. 

In the case of a 100% polarized electron gas, only the 
components of the spin current calculated according to 
Eqs. (33)-(35) corresponding to the majority electrons, 
are relevant. The transition to this case implies the sub- 
stitutions fef — > i/c-f . 

In this section we calculated the wave functions of elec- 
trons in the magnetic wire with a sharp domain wall 
and a localized magnetic moment. The eigenfunctions 
of the Hamiltonian (19) correspond to the spin-polarized 
incoming waves (spin up and down). Obviously, an 
arbitrary-polarized incoming wave is not the eigenfunc- 
tion of the Hamiltonian. Nevertheless, we can still con- 
sider the scattering of electron waves with different spin 
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FIG. 1: Schematic picture of the domain wall. 



polarizations. For example, such a state can be created 
by means of an injection from the tip, and its lifetime 
t can be long enough on the scale of the characteristic 
time of the domain wall motion. In this case, a super- 
position of states with different incoming spin-polarized 
waves can be used to calculate the torque. 



C. Magnetic wire with a thin domain wall 

In the case of a thin metallic wire, when kp-^^d <C 1 
(d is the wire diameter and fepfj are the Fermi momenta 
of minority and majority electrons, respectively), we can 
assume that only one quantization level is filled with elec- 
trons. 

Let us consider a magnetic wire with a single domain 
wall, with the magnetization M being directed along the 
x axis for x < —w and along the opposite direction for 
x > w. Here 2w is the domain wall width, and we chose 
the spin coordinate system like in the previous section, 
with the x axis along the wire. 

Upon applying a small voltage an electric current flows 
in the wire. For definiteness, we assume the current to 
flow in the direction opposite to the x axis direction (i.e., 
the electron flux is oriented along x). If the only imper- 
fection in the wire is the presence of the domain wall, 
one can assume a jump in the potential at the wall, 
and both the charge and the spin currents can be cal- 
culated as integrals over energies in the interval between 
Efr and Sfl — £fr + e A0, where efl and Sfr are the 
Fermi levels on the left and right sides. In the limit of 
small voltage, |e A(f>\ <C sp, the transport is linear and is 
associated with electrons at the Fermi level. 

We assume the electrons approaching the domain wall 
from the left are spin-polarized according to the magneti- 
zation direction in the left part of the wire. The incoming 



electrons are scattered from a large number of magnetic 
moments in the domain wall. We consider this scatter- 
ing using the point interaction of electron with each of 
the localized moments. This corresponds to the picture 
with an array of well-separated magnetic moments like 
in the case of magnetic semiconductor doped with mag- 
netic impurities. We assume that there is a large number 
of magnetic atoms with different orientation of moments 
within the domain wall. Accordingly, the electron trans- 
mitted through the wall is multiply scattered from many 
magnetic moments. 

Thus, to calculate the transmission of electrons 
through the domain wall, we should take the pertur- 
bation created by the total magnetic moment M(ir) = 
^ Mji(i — Xi), where M.; is the localized moment at 
the point x = Xi, and all of the moments are lo- 
cated within a region of the domain wall width, < w, 
which in turn is assumed to be small as compared to the 
wavelength of electrons, kp^^m <C 1. 

The scattering from the total moment 'M.(x) located 
within a region much smaller than the electron wave- 
length can be described using the model of a spin- 
dependent delta-function potential, 11 " Then, in the limit 
of small voltage, we obtain the current 



.in 



; A</> 



2nh 



T/ 



+ l%l 2 + ^. 



(39) 

where the tilde means the transmission coefficients for 
the scattering of electrons from an effective momenliAi 
M e /y ~ f^™ M(x) dx. It is the Biittiker-Landauer for- 
mula for conductivity, which can be obtained in the lin- 
ear response approximation, using the basis of scattering 
states. There are two contributions in Eq. (39) related 
to the incoming waves with spin up and spin down po- 
larizations, and with the corresponding momenta at the 
Fermi surface, fcf^ = fepT,!- 

In the case the domain wall possesses a magnetization 
profile as depicted in Fig. 1, the effective moment M e // 
is oriented along the y axis. The transmission coefficients 

(22), (23) and 
= 1, and the 



t|, i|/ and ti, t^f have the form of Eqs. 



(31), (32), respectively, with n x 



0. 



substitution go — > go = 2mgM e f f /h . We can relate the 
magnitude of M e ff to the continuous magnetic profile 
within the wall, M e ff ~ J_ M y (x) dx. 

The spin current can be also calculated in linear re- 
sponse approximation using the scattering states^i It 
includes the sum of partial spin currents 



e A(j> ( }%{x) i s Ax) 



2wh 



v 1 



"1 



(40) 



where the 



components of j| ^ can be found using 
Eqs. (24)-(26) and (33)-(35)' with the substitution 
*T,ii h,U ^ h,l> h,lf- The appearance of w T and Wj. in 
the denominators of Eq. (40) is related to the ID density 
of states for spin up and down electrons. The compo- 
nents of the spin current perpendicular to the x axis, are 
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oscillating functions. As we see from Eqs. (24)-(26) and 
(33)-(35), the wavelength of oscillations is determined by 
the inverse momentum at the Fermi level. Hence, the 
wavelength of oscillation of the transverse component of 
the spin current is much larger than the domain wall 
width. 

It is worth noting that in three-dimensional systems, 
the transverse component of the spin current decays due 
to the integration over momentum in the DW plane. In 
metallic ferromagnets, the decay is very fast due to the 
large Fermi momentum of electrons. However, there is a 
nonvanishing spin transfer for the transverse component 
in the 3D case, too. 

We can also calculate the accumulated spin density in- 
duced by the external current jo- It can be found as the 
expectation value of the spin <7 M in the scattering state 
of the incoming electrons, integrated over all energies be- 
tween Ep and ep + e Acj>, similar to the calculation of the 
charge and spin currents. Accordingly, we find 
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2irft, \ Vf v± 
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where S|^(0) can be found using Eqs. (36)-(38) with 
n x = and the substitution go — > go, which corresponds 
to the scattering from the effective magnetic moment 

Me//- 

Finally, we find the torque acting on a single localized 
moment in the domain wall. For this purpose we use 
Eq. (16) with S(0) from (41), describing the spin accu- 
mulation created by scattering from the domain wall as a 
whole. In its turn, Sf^O) is calculated as explained after 
Eq. (41) using (36)-(38). The result can be presented in 
the general form 



T = — [n n x (nxs)l(nxsl 

e 



where 



9o 9o (kf 



2fe T A; i (A;t+A; i ) 2 +^(A;f 



goi^ + ki) 2 [(fc T 



g 2 o] 



(42) 



(43) 



(44) 



and s is the unit vector along the spin polarization cor- 
responding to magnetization M at x < —w. The depen- 
dence of the coefficients ry and £ on the parameters go and 
on the electron gas polarization P = — fcj)/ (fef + fc| ) 
is presented in Figs. 2 and 3. As we see, both coefficients 
strongly depend on the parameters of the ferromagnet 
and on the parameters of the wall. In the case of small 
coupling go, we obtain £ 3> r), i.e., the torque is mostly 
related to the second component in Eq. (42). In contrast, 
if go is larger, the first term in (42) dominates. 



FIG. 2: Dependence of the factor r\ on the effective coupling 
go for different values of the electron polarization P. 




FIG. 3: Coefficient £ vs. 
values of P. 



coupling constant go for different 
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D. Spin torque in p-type magnetic semiconductors 

The conductivity of magnetic semiconductors like 
GaMnAs is usually of the p-type. The valence band 
of these semiconductors can be described by a matrix 
Hamiltonian, which includes the spin-orbit interaction^ 
Thus, the calculation of the transmission of holes through 
the domain wall in GaMnAs semiconductors needs a dif- 
ferent model, which takes into account the complex band 
structure. 

Here we use the Luttinger model for the energy spec- 
trum of holes with the momentum J = 3/2^4 We also ne- 
glect the anisotropy of the energy spectrum. To simplify 
the calculations, in this section we take the quantization 
axis z along the wire. Then in the quasi-one-dimensional 
case with the domain wall in the y — z plane, the Hamil- 
tonian of holes is 



H 



2mo 



572 \ d 2 h 72 
7i + — ) — 



J 2 



dz 2 



2 / dz 2 mo 
■g[J y M y (z) + J z M z (z)}, (45) 



where mo is the mass of free electron, 71 and 72 are the 
Luttinger parameters, and J M are the matrices of the to- 
tal momentum 3/2. Note that we are using Hamiltonian 
(45) to describe the holes like unfilled electron states in 
the valence band. Then the correct statistics of holes cor- 
responds to the reversed energy axis of holes as compared 
to that of electrons. 

Similar to the previous consideration, we take the mag- 
netization M along the axis z for z < — w and in the 
opposite direction for z > w, whereas in the region of 
— w < z < w the moment changes its orientation rotating 
in y — z plane. For z < —w the holes can be described by 
the energy spectrum consisting of four parabolic bands 
of the particles labelled by the momentum projection J z 



E 



±3 



/2 (fc) 



h 2 k 2 3gM 
2m, T 2 



E±i/2 (k) 



2 1,2 



2m/ 



gM 



(46) 



(47) 



where m t = m /(7i - 272) and to, = m /(7i + 272) 
are the masses of heavy and light holes, respectively. In 
accordance with Eqs. (46) and (47), the energy band of 
the heavy holes with the moment projection J z = —3/2 
is above all other bands. In the region of z > the 
spectrum is the same but with the opposite signs of J z . 

We assume that the holes are fully polarized so that the 
hole density is rather small. Correspondingly, we assume 
that the chemical potential fj, is located in the interval of 
energies gM/2 < ^ < 3gM/2. 

The scattering state of holes corresponding to the wave 
J z = —3/2 incoming from z = — 00 is 



-ikz , r * e ikz \ 



z < -w, (48) 



^( z ) = (t*e- lkz , t*e- KlZ , t*e' K2Z , 



tie-"**) , 
z > +w. (49) 



where r, r\, ...r^ and t, t\, ...£3 are the reflection and the 
transmission coefficients, respectively. The momentum of 
heavy hole k is taken at the Fermi surface, —h 2 k 2 /2mt + 
3<?AT/2 = s. The other momenta Ki correspond to 
the decaying components of the wave function, k\ — 

[2m, (e - gM/2) /h 2 ] 1/2 , k 2 = [2m, (e + gM/2) /h 2 ] 1,2 , 

1 /2 

and k 3 = [2m t (e + 3gM/2) /h 2 ] . Note that the trans- 
mission coefficient t in this notation corresponds to the 
transmission from the state with moment J z — —3/2 to 
the state J z =3/2. 

In the limit of w — > 0, the matching condition can be 
presented in the matrix form 



diag (m~ t 




dif; 
dz 



-A Jy </>(0) 



-5, 

0, 



(50) 



where Ao = 2gM e ff/H 2 . 

Using Eq. (50) and the continuity of the wave function 
at z = we can calculate eight reflection and transmis- 
sion coefficients. The accumulated spin density S(0) in- 
duced by the current flowing along z can be calculated 
like in the previous section, but with the opposite sign 
because the accumulation of polarized holes means a loss 
of real particles (electrons). Then we find 



S x (0) 



6 -4^- Im (Vstlt + 2«2*i + VzttU 

■KflVt V 



(51) 



S (0) = -^4^" Re (V3*i* + 2i*i! + V3t* 3 t 2 ) , (52) 
2-nnvi V ~ / 



S z (0) 



e 



(3|*| 2 + |ii| 2 -|i2| 2 -3|* 3 



(53) 



where v t = hk/mt is the velocity of heavy holes at the 
Fermi level, e A</> = efr ~ £fl > 0, and efl and £fr 
are the Fermi levels at z < ~w and z > w, respectively. 

Using Eqs. (16) and (42), we find the parameters 77 
and Q determining the torque acting on a single magnetic 
moment Mq 



V 



C 



egMp 
jo ?i 

egM 
jo ft 



S z (0), 



(54) 



(55) 



where jo = — e 2 A</> \ t\ 2 /2nh, and the "-" sign in current 
is due to the opposite charge of holes. 

The dependence of r\ and £ on the magnitude of mag- 
netic splitting gM for different bulk hole densities p is 
presented in Figs. 4 and 5. We take the cross section 
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gM (eV) 



0.8 



FIG. 4: Dependence of the factor r\ on the magnetic split- 
ting gM in the valence band of magnetic semiconductor for 
different values of the bulk hole density p. 



A = 1 nm 2 , and the momentum of heavy holes k = irpm, 
where pm is the linear density of holes. 

As we can see from Figs. 4 and 5, the factor £ is negligi- 
bly small as compared to r\. In our model, the density of 
holes and the spin splitting are independent parameters. 
Thus, the magnitude of torque r\ increases with the de- 
creasing hole density p at a fixed value of gM. However, 
in real magnetic semiconductors these values are not in- 
dependent, and the magnetic splitting increases with the 
increasing hole density^ 

It should be noted that the used condition of w — > 
implies that not only the wavelength of holes with J z = 
3/2 is large as compared to the domain wall width, kw <C 
1, but also the conditions KiW -€i 1 for all Ki should be 
fulfilled. This condition is restrictive for the magnitude 
of the magnetic splitting, {gMmtY^w jfi <C 1. 



III. MOTION OF THE DOMAIN WALL 

A. Hamiltonian and equations of motion 

Now we consider the Hamiltonian H. describing a 
quasi-one-dimensional magnetic system with a domain 
wall. We adopt a model including the magnetic exchange 
interaction and two different anisotropy constants Ai and 
A2 in the z and y directions, respectively (see Fig. 1) 



Hq 



On 

dx 



Ai 2 

2 2 



A 2 2 

Y y 



(56) 
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FIG. 5: Coefficient C, vs magnetic splitting gM for different 
values of p. 



where a is the constant of exchange interaction, and n(x) 
is the unit vector of magnetization. Using this Hamilto- 
nian we are going to describe the magnetic nanowire like 
presented in Fig. 1. Correspondingly, we assume that the 
magnetization vector field n depends only on the coor- 
dinate x and time t. The Hamiltonian TIq is the energy 
density of the magnetic system in the absence of the spin 
torque. 

In the following we restrict ourselves by considering 
the above-calculated spin torque as a driving force act- 
ing on the domain wall. Hence, we neglect the direct 
transfer of momentum from electrons reflected from the 
domain wall. As we show in Appendix, this effect is 
much smaller than the above-calculated spin torque. On 
the other hand, our consideration has a general character 
without specifying the mechanism determining the values 
of the factors ij and 

Using the spherical angles 9{x,t) and ip(x 7 t), we can 
rewrite the Hamiltonian 7io as 



2 \dx 



dtp 

dx 



sin 



Ai 



^2 . 2 
— sin f 



■ cos 



sin 2 <f. 



(57) 



The Landau-Lifshitz-Gilbert equation of motion in- 
cludes a damping term and two possible components of 
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the current-induced torque, as discussed in the previous 
section, 



1 dn 

7 ~dt 



SHo d SHo 



On 



Sn dxS(dn/dx)J dt 

+ JoC n x s + JqT] n x (n x s) , (58) 



where a is the damping constant, 7 = gfiB/hM is the gy- 
romagnetic ratio divided by M, Jo = joh/egflo, and Qq 
is a volume per magnetic moment. In Eq. (58) the spin 
torque is expressed in terms of transferred moment per 
unit volume, and it enters directly the equation of mo- 
tion. The corresponding spin-torque terms in the mag- 
netic Hamiltonian can be presented as 



T-tint = JoC n ■ s + J r) I dm- ( ^ x s 



(i 



(59) 



where n(r = 0) = and n(r = 1) = n. 

The Lagrangian of the magnetic system contains a 
term with a time derivative as follows 9 ' 20 



£ = A dx 



I^(cos0-l)-tt 
7 dt 



(60) 



Neglecting the damping term, the Landau-Lifshitz equa- 
tions of the magnetic dynamics take the form 

1 06 d 2 tp . „ A . „ . 

= —a — — — sm9 + A2 sm9 simp cos tp 



7 dt 



sin 9 dtp 



dx 2 



-JoT] cos 9 cos ip, —JoC sin tp, (61) 



7 



d 2 9 

dx 2 



dip 

dx 



smtf cost 



+Ai cos 9 sin 9 — X 2 sin9 cos 6* sin 2 ip 
+Jorj sin tp + JoC cos 9 cos tp. 



(62) 



In the absence of current, jo = 0, they have 

the well-known^ kink-like static solution ipo(x) = 

1/2 

arccos [tanh(/?oa;)] and #0 = 7r /2, where @o = {M/o) 
is the inverse width of the static domain wall. From now 
on we assume for dehniteness that Ai > A2, so that the 
static domain wall with the magnetization in the x — y 
plane is energetically more favorable. 

In a general case, the solution of the nonlinear dynam- 
ical Eqs. (61)-(62), describing the moving domain wall, is 
a difficult problem. Therefore, we assume in the following 
that one of the anisotropy constants is large, Ai S> A2. 



write 9(x,t) — ir/2 + xi x ^) and take \x{x, t)\ <C 1. Up 
to the second order in x(x,t) field, the Lagrangian C is 



£ = A dx 



7 at 2 \dx J 2 \dx 



(1 2\ Ai 2 ^2.2/1 2 

(1 - x ) - y x - y v ' ~ x 



+ JoV Xsin<£ 
-JoC cos<^] , (63) 



Since we restricted the treatment to quadratic terms 
in x in the Lagrangian, the integral over x is Gaussian, 
and we can integrate ouli^ the x fields to obtain 



£ 



A dx - 
J [2 

1 d(p(x) 



7 



dt 



dx' G(x,x') 
JoT] smip(x) 



1 d(p(x) 
1 



dt 



dtp 
dx 



JoT] sm.{p{x) 

2 \ 

^2 . 2 

sin tp 

2 v 

-JoC cos <p], (64) 



where the Green function G(x, x) obeys the following 
equation 



d 2 ( dtp , . ., 

- a d^- a Yd-x I + A *- A2Sin * 



G[x,x') 
= 5{x-x'). (65) 



Note that ^-fields are taken at the same time t in 
Eq. (64). It follows from the equation for Green func- 
tion G(x, t; x',t') ~ S(t — t') describing the propagation 
in time. 

Equation (64) contains ip(x,t), which should be the 
saddle point solution of the Lagrangian, i.e., the self- 
consistency should be preserved. 

Neglecting the first term in Eq. (65) we can find an 
approximate formula for the Green function proportional 
to 5(x — x') 



G(x, x') = 5(x — x) 



dp 
dx 



+ Ai — A2 sin tp 



(66) 

This form of G(x,x') leads to the point interaction of 
the tp- fields in the first term of Eq. (64). Physically, ne- 
glecting the first term with derivatives in Eq. (65), we 
substitute the finite-range interaction by the 5-like one. 

One can estimate at which conditions the use the 
Green function in the form of Eq. (66) is justified. The 
exact solution of Eq. (65) can be presented as 



G(x,x')=J2 



£n + Al 



(67) 



B. Strong easy-plane anisotropy 

We consider the case of a large easy-plane anisotropy, 
and, accordingly, assume that for the moving domain wall 
(subjected to the torque) the deviation of magnetization 
vector M from the x — y plane is small. Then we can 



where <j) n (x) and e n are the eigenfunctions and corre- 
sponding eigenvalues of the equation 



d 2 ( Otp\ . , 

~ a dx^~ a > ~ A2Sm V 



4> n (x) = 0. 



(68) 
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We expect that the function ip(x) in Eqs. (65) and (68) 
is similar to the form of the static solution ip (x). Thus, 
Eq. (68) corresponds to the Schrodinger equation for a 
particle of mass m = ft 2 /2a in the potential well V(x) 
of width Lq ~ (a/^) 1 / 2 - The energy spectrum of this 
problem consists of a level in the well eo — — A2 and a 
continuous spectrum for all positive energies. 

For (p(x) = ipo(x), the potential has the form of V(x) = 
—2X 2 / cosh 2 (/3 a;), and the discrete energy spectrum^ 
has one level e = — 4A 2 . The eigenfunctions <j) n {x) cor- 
responding to the continuous spectrum, are oscillatory 
functions, so that their contribution to Eq. (67) can be 
estimated as G^ cont \x, x') ~ (aAi)~ 1 / 2 e- K1 ^"^'l , where 
«i = (Ai/a) 1 / 2 . Thanks to Ai ^> A2, it is a strongly 
localized function on the scale of the distance (static 
domain wall width) Lq. On the other hand, the con- 
tribution of the localized state gives us G^°'(x,x') ~ 
(I/L0A1) e~ K °\ x ~ x I, where Ko = l/^o- Thus, in the case 
of Ai ^ A2 (i.e., strong in-plane anisotropy), the con- 
tribution of G<°) can be neglected as compared to the 
short-range interaction. Using the condition of strong 
easy-plane anisotropy, Ai 3> A 2 , we can simplify Eq. (66) 
essentially and obtain the expression 



G{x,x') 



5(x - x') 
Ai 



(69) 



In this approximation, the Lagrangian (64) acquires the 
following form 



C = A 







j dx 


2Xi ( 





1 dp T . A 2 a ( dp x ' 

-—-Jon an?) 



^2.2 1 j j- 

-— sin ip + JqC cos(p 



(70) 
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FIG. 6: Dependence of the parameter f3 on the domain wall 
velocity v for different values of the current. 



The dependence f3{v) is presented in Fig. 6 where we 
denoted v — v/^\fX\a and jo = Jari/\JXiX 2 . 

When jo = 0, we find from (72) that (3 2 = /3 2 /(l - 
v 2 /j 2 Xia). It means, that in the absence of the current, 
the solution for a moving domain wall is more sharp as 
compared to the static wall. 

For jo 7^ and v — > 0, the value of (3 depends on the 
current as f3 2 = 0q (1 — Jfirfa/ XiX 2 ), i.e., the current 
makes the thickness of the static domain wall larger. 

Now we can use the velocity v as a variational param- 
eter to minimize the Lagrangian 



C = ^y I ™r ,-(.<■) ,/(.<, ) 



(73) 



The corresponding saddle-point equation is 



d 2 ip J$r) 



2.2 Qj2 

7 2 Ai BP • X, ^™V + a g^ 
— A2 sin ip cos ip — JqC sin tp = 0. 



(71) 



Now we can study the problem of the domain wall dy- 
namics in terms of a single <p(x, t) field. 



C. Solution for ( = 

Let us consider the possibility of the kink-like solution 
moving with an arbitrary constant velocity v, p(x,t) = 
tp(x-vt). We can find such solutions in the case of ( = 0, 
trying a function which obeys the equality dp{x)/dx = 
f3 sinp(x). This function differs from the static solution 
only by a different choice of (3 instead of /3 = {X 2 /a) 1 / 2 . 
Substituting it into (71), we obtain the equation that 
relates the values of /? and v as 



ii 2 



7 2 Ai 



-A, 



72 2 



= 0. 



(72) 



where 



F(v) 



1 

a7 



v(3{v) 



a/3 2 -X 2 , (74) 



the function (3(v) is defined by Eq. (72), and the integral 
in (73) does not depend on v. 

Using Eqs. (72) to (74) we find that for j = 0, the 
quantity F(v) — — A 2 for any v. Thus, the minimum 
of C corresponds to (3 = (3q, which is the minimum 
value of the dependence (3(v) for jo = 0. In the limit 
of a small velocity, v 2 <C 7 2 Aia, and using the rela- 
tion J sin 2 p(a) da = 2, we find the kinetic energy of 
the moving domain wall in the form of Eki n = m*v 2 /2, 
where m* — A^/X^/^Xi^/a is the effective mass of the 
domain wall. This is in agreement with the definition 
from Ref. for A : = 2ttM 2 . 

For jo ^ 0, we can present the dependence of the factor 
F on both parameters v and jo as 



F(v) = -X 2 



l-2j 2 -2j £ 



1 



~2 \ V2- 

Jo 



(75) 
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In the limit of v 
jo > iocr) where 



0, the factor F changes its sign for 



Jo < 



V2hr] 



(76) 



is the critical current. Thus, if jo > jocr, the solution 
with moving domain wall is energetically favorable. 

We can interpret the effect of the current as leading to 
an effective reduction of the effective mass of the domain 
wall. For jq > jo cr the current induces an instability 
towards a spontaneous motion of the wall. 



D. Case of < / 

In the case of ( ^ 0, there are no solutions of Eq. (71) 
corresponding to the motion of the domain wall with a 
constant velocity because the last term in this equation 
acts as a force accelerating the domain wall. Indeed, 
if we assume a probe solution in the form of <p(x, t) = 
ip [x — xo(t)] we find 



m* x (t) + J A( = 0, 



(77) 



where m is a constant in the limit of a small velocity. 
In other words, Eq. (77) describes the acceleration of the 
domain wall just after we apply some voltage. Hence, our 
model can describe the steady state if we include a viscos- 
ity (friction) into the equation of motion. We can use the 
damping term from Eq. (58). Writing the corresponding 
additional term in Eq. (71) as Fd = —adtp/dt, we find 
the following equation that determines the velocity of the 
moving domain wall as follows 



vf3{v) 



Jot 
a 



(78) 



This equation indicates a linear dependence of the veloc- 
ity on the current in the limit of a small velocity, when 
P is constant. As we see from Eq. (78), it corresponds to 
a large damping. 

Effective friction may as well stem from the pinning 
by impurities. This case can be described phenomcno- 
logically leading to another mechanism for the critical 
current & 



IV. CONCLUSIONS 

We calculated the components of the spin torque acting 
on a thin domain wall in a magnetic nanowire subject 
to an electric current. These components can induce a 
rotation of magnetic moments in different directions. 

We also considered the dynamics of a domain wall in 
the presence of the current. We demonstrated that a 
moving magnetic kink, similar to the static domain wall, 
can be a solution of the equations for the magnetic dy- 
namics only at some special conditions. We found these 



conditions in the case of a large ratio of the magnetic 
anisotropy constants. In the limit of small velocities, the 
solution does look like a kink but its width decreases with 
increasing velocity. In the limit of a small velocity, the 
domain wall moves like a particle of a mass determined 
by the exchange interaction and anisotropies. One of 
the spin torque components dominating at the small 
coupling, acts as a moving force on the domain wall, ac- 
celerating the wall, provided that there is no pinning to 
impurities. 

Recent direct observations of the domain- wall configu- 
rations show that the spin structure of the wall changes 
with the current, and the structure depends on the ve- 
locity of the domain wall motion^ 

We performed the calculation of torque in the limit of 
thin domain wall, w <C \f- ft allows to simplify essen- 
tially the problem and to do all the calculations analyti- 
cally. In reality, this inequality may be not well satisfied 
even in the magnetic semiconductors. Here we present 
the estimations for the case when the above-mentioned 
inequality is obeyed. 

Let us take the cross section of the wire A = 1 nm 2 and 
the bulk carrier density n^o = f0 19 cm~ 3 , corresponding 
to the linear density n\u — n^jjA — 10 5 cm -1 . It gives 
us kp = tttiid ~ 3 x 10 5 cm -1 , and the carrier wavelength 
\ F = 2Tr/k F ~ 100 nm. 

To estimate the domain wall width, we assume the 
magnitude of magnetization M — 100 Oe, the demagne- 
tizing factor along the y axis — 0.3, and calculate the 
anisotropy constant as A2 ~ 87m M 2 ~ 10 5 erg/cm 3 . 
For the energy of magnetic interaction Ei nt ~ 10 meV 
at a distance between magnetic ions of cq = 1 nm, 
the exchange parameter a — E; mt cq/A ~ 10~ 8 erg/cm. 
Then the domain wall width has a reasonable value of 
w = (C1/X2) 1 ' 2 — 10 nm. Comparing these estimations, 
we see that the main inequality of w <C Xf is satisfied. 
At a larger carrier density, both w and A can be of the 
same order of magnitude or the inequality is reversed like 
in magnetic metals. In this case, the constants £ and r\ 
should be calculated numerically. 
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APPENDIX A: SPIN TORQUE DUE TO THE 
MOMENTUM TRANSFER 

The reflection of electrons from the domain wall is ac- 
companied by the transfer of momentum from the elec- 
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tron system to the domain wall. In the presence of elec- 
tric current transmitted through the magnetic wire, it 
creates an additional force acting on the walli& Here we 
estimate the magnitude of this effect in the case of a thin 
DW, k F w < 1. 

The force F is determined by the total transferred mo- 
mentum in unit time. Taking into account the contribu- 
tion of spin up and down scattering states corresponding 
to the waves incoming from — oo in the energy range be- 
tween e p and e + e A.<fi, we find 



T - 



eAcj) 
~2^T 



1+ 1 2 I |2 V] ,, ,2 

l*T/l +-I r i/I 



+fc i (l + K| 2 -|t i/ | 2 + ^|r T/ | 2 -^|t T | 2 



.(Al) 



This force tends to shift the domain wall along the x 
direction. For a local moment within the wall it is equiv- 
alent to the appearance of a torque. To estimate the 
magnitude of this mechanical torque acting on a single 
moment we use a simplified model. 

We describe the domain wall by the ip(x) field, which 
is the angle in x — y plane determining the direction of 
moment M(a;) as shown in Fig. 1. We assume that the 
shift along x is related to the following interaction 



Ti-int = X(p(x) v(x) 



(A2) 



where A is a constant, v(x) = —depo/dx, and <fo(x) is 
the static solution for the domain wall. The potential 



v(x) has the form of a potential well in the vicinity of 
the domain wall, and it forces (makes energetically fa- 
vorable) a correction to the <f{x) field of the same form, 
Stp(x) ~ dipo/dx. On the other hand, the correction 
5ip(x) = {dip/ dx) 5xq is the shift along the axis x by 
8xq. Thus, the interaction term in form of (A2) in the 
equation of motion for the <p(x) acts as a shifting force. 

The constant A should be determined by the condition 
that the energy 5E associated with the shift, gives the 
force T: 



T 



SE 
Sx 



XA 



V dx 



dx. 



(A3) 



Using the known solution, dipo/dx = /3ainipo(x), we find 
X = F/2{3A. 

The equation of motion for (p(x) (64) includes the ad- 
ditional torque term as Xv(x) — Tv(x) /2/3A. Using (Al) 
we estimate the torque acting on the localized moment 
M = MQ 



e A ' 



(A4) 



where is the volume of an elementary cell. We find 
that the relative contribution of the momentum-induced 
torque with respect to the spin transfer is 



T mt /T st ~ k F Q/A < 1. 



(A5) 
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